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Abstract 

The viscous transonic flow around a low- aspect-ratio wing has been computed 
using an implicit, three-dimensional, “thin-layer” Navier-Stokes solver. The grid 
around the geometry of interest is obtained numerically as a solution to a Dirichlet 
problem for the cube. The geometry chosen for this study is a low-aspect-ratio wing 
with large sweep, twist, taper, and camber. The topology chosen to wrap the mesh 
around the wing with good tip resolution is a C-0 type mesh. Using this grid, the 
flow around the wing was computed for a free-stream Mach number of 0.82 at an 
angle of attack of 5°. At this Mach number, an oblique shock forms on the upper 
surface of the wing, and a tip vortex and three-dimensional flow separation off the 
wing surface are observed. Particle path lines indicate that the three-dimensional 
flow separation on the wing surface is part of the roots of the tip-vortex formation. 
Tne iming or tne tip vortex berore tne wing trailing edge is clearly observed by 
following the trajectory of particles released around the wing tip. 


I. Introduction 

It is well known from classic wing theory 1 that the presence of a vortex wake 
is a necessary condition for all finite-span lifting wings. The most advanced full- 
potential numerical methods for the simulation of subsonic and transonic flows com- 
bine classic wing theory (vortex sheet trailing the wing at a fixed spatial location) 
with the full-potential calculation to compute realistic lift values for wings of fairly 
general shapes. 2 In reality, this vortex sheet will roll up into discrete vortices. 1 This 
roll up starts close to the tip of the wing where the flow is predominantly three- 
dimensional and possesses a concentrated (rolled-up) vortex structure. The wake 
vortex generated by large transport aircraft is known for its potential hazard to 
smaller aircraft. This hazard is the major impediment to increasing traffic capacity 
at major air terminals and a large effort is being directed to minimize it. 3 In the 
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case of low-aspect-ratio wings, the tip vortex contributes directly to lift and drag 
and, hence, it is important in practical applications. 

To contribute to our understanding of the formation of the tip vortex, a 
mathematical model is needed in which the appearance of the tip vortex is a 
result of the solution rather than an empirical input to the solution. At high 
Reynolds number, the simplest model seems to be the thin-layer approximation 4 
to the Navier-Stokes equations, since the vortex formation is a viscous phenomenon 
generated in a thin layer close to the surface of the wing, whereas most of the flow 
field is convection dominated. The model that incorporates the full viscous terms 
is preferable when a fine mesh is used in all three-directions. However, experience 
with two-dimensional computations indicates that the difference between the results 
from a full Navier-Stokes calculation and a thin-layer calculation is negligible when 
a coarse grid is used in the streamwise direction. The objective of this paper is to 
study the formation of the tip vortex off a low-aspect-ratio wing. The thin-layer 
approximation is used as a first model for this purpose. Vortex formation and 
capturing at wing tips are reported by Rizzi and Eriksson 5 and Eriksson and Rizzi 6 
who use the Euler equations as the mathematical model. It is not clear why a tip 
vortex is observed in inviscid (Euler) flow computations, and the understanding of 
this observation is still under investigation by Rizzi. 7 

Considerable progress has been made in developing the solution methods for the 
nonlinear equations of fluid motion. For a review of the achievements of the numeri- 
cal methods in computing transonic flows, the reader is referred to the recent article 
of Mehta and Lomax. 8 The methods range from the fully explicit to the factored 
fully implicit schemes. The application of these methods to three-dimensional flows 
has been restricted mostly to simple geometries. Both the computational efforts 
involved and the problem of grid generation over complicated geometries seem to 
hinder the application of the available computational methods to three dimensions. 

The procedure followed in this work is to treat the mesh-generation problem 
and the solution to the governing equations as two separate operations. This appears 
to be the approach followed by most authors 6 ’ 0,10 to compute the three-dimensional 
flow over complicated geometries. The adaptive gridding approach 11 is still under 
development, and its application to three-dimensions for complicated geometries is 
still not proved. The mesh-generation method chosen to generate the grid around 
the geometry of interest is the elliptic mesh generator described in Sec. H. The 
geometry chosen for this study is a low-aspect-ratio (A r — 0.85) wing with 45° 
sweep, 8.3° twist, and 0.3 taper. The mathematical model used to compute the flow 
around this wing is the thin-layer Navier-Stokes equations in generalized coordinates 
(described in Sec. III). The numerical method used to solve these equations is 
the factored implicit scheme 12,13 as applied to three-dimensional aerodynamics 
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applications by Pulliam and Steger. 14 The data-management strategy described in 
Sec. IV allows the use of fairly fine three-dimensional meshs (up to 200,000 grid 
nodes). The computed results (using 75,000 grid points) for the formation of the 
tip vortex over a low-aspect-ratio wing at a transonic Mach number (M = 0.82) 
are presented in Sec. V. 


n. Grid Generation 

The grid around the geometry of interest is obtained numerically as a solution 
to a Dirichlet problem for the cube. The approach relies on the method of Mastin 
and Thompson, 15 with the source terms specified using a variation of the method 
of Thomas. 16 The source terms in this approach are specified so that the Poisson 
equations possess exponential solutions for locally constant parameters. This allows 
grid control inside the volume by the grid-point distribution assigned on the surface 
of the domain. The distribution on the surface of the domain is controlled by 
the grid-point distribution along the boundary of the surface using the method of 
Thomas and Middlecoff. 17 

Generation Procedure 

The first task is to distribute the grid points on the boundary of the six surfaces 
of the volume of interest. This is done using a composite transformation consisting 
of a hyperbolic-tangent and a hyperbolic- sine functions as follows (T. Morel, private 
communication, Integral Technologies, Westmont, 111., 1983.): 


a{x — Xq) = sinh( cT + d) (1) 

ee + / = tanh~ 1 ( ! ,r + />) (2) 

where x is the physical coordinate € [x\,X 2 ], £ is the transformed coordinate 6 
[1, JV], and T is an intermediate coordinate 6 [1, AT]; N is the total number of grid 
points. To determine the seven constants a, c, d, e, f , g, and h, the following four 
conditions are imposed: 

Y— i} w ^ en x ~ Xl > ancl pZ jv} w h en x — x % ( 3 ) 

We are left with three adjustable stretching parameters that we specify as follows: 

h — 9 + h 

te = gN + h (4) 

t 3 — cN + d 
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Parameters t\ and t 2 provide grid clustering around x\ and X 2 , and t% provides grid 
clustering around so- 

Once the distribution of the grid points on the boundary of the surface has 
been determined, a two-dimensional elliptic grid generator is used to distribute the 
grid point on the surface. The generating system is of the form 

<*(*« + <pxz) - 2pXSt, + l{X V r, + tXq) = 0 

a{y^ -f- 02/?) — 20Vtn + 7 {Vm + i’Vv) = 0 1 ^ 

where 

a = x* -f y v 

P — XfXt, + y^y„ (6) 

7 = *1 + 2/? 

Imposing the orthogonality condition at the boundary on the grid lines yields the 
values of 0 and ip on the boundary as follows: 17 

0= + ftS/«)/(®? + 2$ on r) = l,Nr, . 

0= —{XtjXnn + y v yr)r>)/( x v + y 2 v ) on Z — l > N Z 

The values of the parameter 0 (or ip) at the interior are then computed by linear 
interpolation along £ =const ( tj =const) lines. 

In the above formulations, the boundary is assumed to be rectilinear and the 
surface is assumed to be planar. For general lines (or surfaces) the grid is generated 
on an appropriate rectilinear (or planar) projection, and the points on the line (or 
surface) are deduced from the points on the projection. A more general method 
that accounts for the effects of curvature would require specifying the tangents and 
radii of curvature, which in general may not be known analytically. 

Once the grid is generated on the six surfaces of the volume, the mapping 
between the physical space (x,y, z) and the computational space (£,??, f) is deter- 
mined by solving the following set of transformed Poisson equations: 

+ 4> x d + {Ztw -j- ipx t ,) + a3(a: ff + ux q ) + 

2{Pi x Sr) + P2 x ti$ + P& x st) — 0 

a l(2/?? 4 ~ 02/?) + a 2(yr)Ti + 02/17) + »3(2/ff H~ ^ y f ) + /o\ 

2(/?i2/?*? 4* feytis 4~ /?32/j?) — 0 

<11(2^ + <pZf) + <*2(21717 + 0^77) + Oi 3(2^ + + 

2(01 ~f~ fozr/f + PaZrf) = 0 
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where 


<*i “ S* H" (y + £« 

<*2 = nl + Vy + Vz 

01 3 — ?* .+ <Ty 4- d 

/?1 = 6cfyr + tyVy + 

@2 == »?3:Ca; + tyyfy + *?*£* 

Pz = &€s + fy£y + UU 

The metric terms £x,£yr • • are expressed in terms of x^,x Vr . . as follows: 

& = J{y n *s - y< z v) v* = J{z & ? - ^) 

fy — — ^ij) *?y = J{XfZ t - X f Zz) 

Zm = J(x n y s - x^) ri z = J(ysx t — y { xf) 

Cx = J(y^ v - ^) 

fy = JU^Xr, — 2 ^) 

= Jixtft, - x n y^) 


( 11 ) 


and 

J~ l = x^y^ + x^Zf, + x n y^ 

x iVv z i~ x vV^ 

Direct extension of the two-dimensional formulation [Eq. (4)] to three dimensions 
yields the values of <j>, ip, and co on the boundary as follows: 

0 = —(£{£(£ + yzyiz + z^)/(x | + y\ + z\) 

on i) = 1, N v and $ — 1, iVj 


0 = ~{x v Xr, v + yr>yrin + z n z r,v)/(tf H- y% + *?) 

on f = 1, iV$ and f — 1, N f 


CO 


= -(z { x s{ + y<y<s + z<z«)/(x* + 2/? H- **) 
on £ = 1, iV^ and r/ — 1,N^ 


( 12 ) 


The continuous representation of the parameters is found in planes f =const, 
T) =const, and f —const by solving analytically the following equations: 


— 0 

(13) 

o 

The solutions to Eqs. (13) are elementary extensions to two-dimensions of the one- 
dimensioDal linear interpolation method. 
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Grid Type 


The geometry chosen for this study is the low-aspect-ratio wing (wing C) 
investigated experimentally and numerically by Hinson and Burdges. 18 The same 
geometry was also investigated experimentally by Keener (private communication, 
NASA Ames Research Center, Moffett Field, Calif., 1982.) and numerically by 
Holst et a/. 19 The characteristics of the wing are described in Fig. 1: 0.85 aspect 
ratio, 45° sweep, 8.3° twist, 0.3 taper. 

In these numerical studies, the full-potential equation was used to compute the 
flow around the wing. With the full-potential methods, the vortex sheets are fitted 
rather than captured, and, therefore, the wrap around the wing tip is not critical. 
If the flow around the wing tip is of interest, then a mesh that wraps around the 
tip is desirable. Different grid wraps can be used to resolve the flow field around 
the wing that will achieve for a given total number of grid points concentration of 
grid points in different regions of the field. Rizzi and Eriksson 5 used the density 
of grid points around the wing as a criterion by which to judge the different grid 
wraps and concluded that the 0-0 and C-0 are the two “best” types. The C-0 
type mesh shown in Fig. 2 was chosen in this study. 

As pointed out by Rizzi and Eriksson, in order to map the wrap into a box-like 
region we need to introduce artificial singular lines. The location of these lines 
should be away from areas of interest since the flow field in the neighborhood of 
these lines may be in doubt. The map into the computational box chosen for this 
study is shown in Fig. 3. The cuts were chosen to be on the leading edge and 
trailing edge leaving the area around the tip free from these cuts. The wing and 
wake surface is mapped into the “bottom” surface of the parametric box (plane 
f = 1); the outer boundary (except the downstream surface) is mapped into the 
“top” surface (plane f == iV f ); the upper plane of symmetry is mapped into the 
“front” surface of the box (plane t] — 1); the lower plane of symmetry is mapped 
into the “back” surface of the box (plane r\ — N v ); the cut at the leading edge 
is mapped into the “left” surface of the box (plane £ — 1); and the downstream 
surface of the outer boundary plane is mapped into the “right” surface of the box 
(plane f = N(). 


HI. Governing Equations 

At high Reynolds numbers, the flow around the wing can be considered to be 
mostly inviscid, except in a thin-layer region close to the body where the viscous 
effects become important. With this in mind, the thin-layer approximation is used 
to simplify the Navier-Stokes equations in generalized coordinates to yield: 14 
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where 




- P - 



pU I 



pu 



pull 4- i x p 


pv 

f 

E = J~ 1 

PVU 4" ZyP 
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■ e - 



( e + p)u J 


r 
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pwW 4- UP 


(e 4 -P)V J 


. {e + p)W . 

r - 1 X 




0 



+ f J + fJK + + fyVf + &«>*)& 

/^ e (<r| 4- <r 2 4- fJK 4- (a*«/3 x?*« s + sw 4 - &«»*)& 
M«(rl + f J 4- dH + (/<«/3X?*ttf + Cyv s + 

{(d + fj + ^)[0.5p e (u 2 + w 2 4- W 2 ) ? 

+/cPr _ 1 ( 7 - l)*” 1 ^ 2 )*] 

. +W 3)(<r a w + f,v + + fyWj + &«/,)}. 


U — £ X U + + £,w 

V = -f TjyV + r\ z M) 

W = fctt + fct/ + frttf 


(15) 


(16) 


(17) 


where U, V, and W are contravariant velocities written without metric normaliza- 
tion. The above set of equations is nondimensionalized using (the free-stream 
sound speed) for velocity, the cord length (c) at the root of the wing for distances, 
and poo the free-stream density for density. One more equation of state is needed 
to close the above system, and we use the perfect gas law to relate the pressure, 
density, and velocity as follows: 


p = (--y — i)[e — 0.5p(u 2 v 2 + w 2 )] 


(18) 


The effective viscosity in Eq. (16) is set equal to 


p e = 1 + Ht (19) 

where pt is a turbulent viscosity computed using the two-layer algebraic viscosity 
model of Baldwin and Lomax. 4 



IV. Numerical Method 


The approach followed in this work is the same as the one developed by Pulliam 
and Steger. 14 The numerical algorithm is the implicit noniterative approximate 
factorization scheme of Beam and Warming. 12 Fourth-order dissipation terms are 
added to the equations to control nonlinear instabilities, and implicit second differ- 
ence terms operating on Aq are added to extend the linear stability bound of the 
fourth-order terms. 

Linearization and Approximate-Factorization 

The equations of motion are differenced using a standard (non-staggered) mesh, 
central differencing in space, and backward (fully implicit) differencing in time. The 
differenced equations are linearized in time, then factored using the Douglas and 
Gunn 20 technique to reduce the three-dimensional problem to three one-dimensional 
problems at a given time-level. The differenced equations are written as follows: 


D^DfjD^Aq = H 

where 

Aq= q n + x - q n 

Dz= r+h6z(dE n /dq n )- erJ-'Vf A^J 
D v = I + h6 r) {dF n /dq n ) - erJ- 1 Vr,A v J 
D f = / - f hdf(dG n /dq n ) - hRe~ 1 8^dS n /dq n ) 

-e/Z-^Vj A f 7 

H = - At{6^E n + 8 v F n + 6 S G n - Re-^S”) 

-es.l-'KV^) 2 + (V,A,) 2 + (V f A f f]Jq n 

The 0 ’s are central-difference operators, V and A are backward and forward difference 
operators, and h = At is the time-step. For a detailed description of the Jacobian 
matrices dE n /dq n , dF n /dq n , dG n /dq n , and dS n /dq n the interested reader is 
referred to the paper of Pulliam and Steger. 14 

Data Management 

Proper resolution of the flow field around the wing requires a large number of 
grid points. Chapman 21 estimated the following requirements to resolve the three-? 
dimensional boundary layer on a wing: 


( 20 ) 


( 21 ) 


7V C = 4.5n^Re°- 2 

N v = 4.5A r n^ReJ- 2 (22) 

AT f =n f 

where n^, n n , and n f are the average number of grid points per length 8 (the 
boundary-layer thickness) in the streamwise, spanwise, and normal direction, respec- 
tively; A r is the aspect ratio; and Re c is the Reynolds number based on the cord 
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length. For Re c = 1. X 10 6 , A r = 0.85, = 0.4, n n — 0.4, and = 15, we 

need — 30, N n — 30, and N s = 15 to resolve the boundary layer on the upper 
surface. (Chapman recommends using = 1, n v — 0.5, and = 20.) Adding 10 
grid points in the ^-direction to capture the wake, doubling the grid points in the 
^-direction to resolve the lower surface, and adding 15 grid points in the f -direction 
for the outer inviscid field yields a required total number of grid points equal to 
72,000. We need 19 words of memory at each grid point to store the independent 
variables and the metrics; the required total number of words of memory is then 
1.4 x 10 6 . The Crayl-S machine available at Ames Research Center is nominally 
a 1 x 10 6 words-in-core-memory machine. We then need to use the disk to store 
the flow-field data. To carry out the computation, part of the data are transferred 
into the core memory of the machine. The data-management algorithm devised is 
a five-plane algorithm. 

The flow field is stored on the solid-state disk (SSD) and the “in-core” working 
set is a slab of five planes. With this scheme, data are transferred to the core 
of the machine in £-;? planes. When the data are needed in the ^-direction in 
more than five planes (for example for matrix inversion), a sweep through the data 
is carried out. The sweep through the data occurs in the f-direction only. For 
a matrix inversion in the ^-direction, two data-sweeps are needed. The solution 
proceeds as follows: during the first data-sweep, two matrices are inverted in each 
£-?? plane. Also, forward elimination for the third matrix inversion is completed. 
During the second sweep* the inversion of the third matrix is finished and the 
solution is advanced. A mesh up to 200,000 grid points can be handled by this 
algorithm. 

Boundary Conditions 

The physical boundary conditions are the conditions at infinity and the no-slip 
conditions on the surface of the wing. However, additional boundary conditions are 
needed when the equations are discretized. In general, boundary conditions on the 
six surfaces of the computational box are required to close the system of differenced 
equations. The use of appropriate boundary conditions is dependent on the mesh 
wrap. The type of boundary conditions used with the C-0 wrap in our case is 
shown in Fig. 3. At f = 1, “cut” conditions are imposed. The cut boundary 
condition used in the present scheme is that the values at the cut are given as the 
average value of the two neighboring points in the computational domain. At £ = 
N$, exit boundary conditions are imposed where the values at the boundary are 
extrapolated from the interior of the domain along £-lines. At rj = 1 and r) — N v , 
symmetry boundary conditions are imposed. At £ = 1, no-slip boundary conditions 
are imposed on the surface of the wing and cut condit ions are imposed in the wake. 
At f = Nf, inflow conditions are imposed where the velocity, density, and Mach 
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number are specified. 


V. Results 

Using the grid described in Sec. II, the flow around the wing was computed 
for free-stream Mach number M — 0.82 at an angle of attack a = 5°. The design 
conditions 18 for this wing were M — 0.85 and a — 5°. The results presented in 
this section are, hence, for a Mach number slightly below the design Mach number, 
and the flow around the wing is expected to be smooth. (Smooth flow is loosely 
defined in this paper as a flow that does not exhibit a massive separation.) 

The surface-pressure contours (Fig. 4) and the Mach contours (Fig. 5) on a 
surface parallel to the surface of the wing show the presence of an oblique shock 
between y ~ 0.4 c and y ~ 0.75 c at about 0.25 the local cord length from the 
leading edge. This shock is also detected by the Cp profiles of the experimental 
data. In the region of the shock, our mesh is relatively coarse, as can be seen in Fig. 
6 where the open symbols represent the computed Cp results at the mesh locations. 
Comparison of the computed Cp profile with the experimental data shows that 
the suction peak is overpredicted. No Mach-mmiber or angle-of-attack corrections 
(for the presence of the walls in the experiment) were attempted to match the 
experimental data. Using the full-potential-equation formulation, Subramanian et 
a/. 2 found that by correcting the Mach number and the angle of attack they can 
match the experimental data at the different spanwise locations. In general, the 
shape of the Cp profiles predicted by the present computation is similar to the 
experimental profiles at the various spanwise locations, but the suction peak is 
overpredicted. 

The flow pattern over the surface of the wing was visualized experimentally 
by Keener using an oil-flow technique. It was found that for low Mach numbers 
(M = 0.7) the flow is smooth over the wing, but at high Mach numbers (M > 0.85) 
a massive three-dimensional separation is observed on the upper surface of the wing. 
The pattern indicated by the oil flow for M = 0.82 is sketched in Fig. 7 from a 
photograph. The major feature of the flow at this Mach number is the appearance 
of a ridge where the oil path lines change direction. The flow-field pattern from the 
numerical solution can be best visualized by plotting the path of particles whose 
motion is defined by the following equation: 



The above equation can also be used as the definition of the streamlines. If we 
restrict the path of the particles to a surface, we will be defining the streamlines 
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of the projection of the velocity vector on that surface. On the surface of the wing 
the velocity is zero because of the no- slip condition and, hence, particles on the 
wing surface do not move. At one mesh ( f = 2 ) above the surface, the velocity 
vector is nearly parallel to the surface of the wing, and by plotting the particle 
path in that plane we will be observing the streamlines close to the surface, but 
not simulating the oil-flow pattern. In general, the oil film has a thickness, mass, 
and surface tension, and a complete simulation of the film is much more involved 
than massless particle tracing. Figure 8 shows the particle path pattern in the 
f = 2 plane (average distance to the wing d 2 — 0.4 X 10“ 3 c). The appearance 
of a ridge is clearly seen where the paths of the particles converge. Converging 
lines indicate flow separation, and diverging lines indicate flow reattachment. Both 
separation and reattachment can be detected in Fig 8. If we compare the computed 
streamlines in that plane (f = 2) with the experimental oil-film pattern, we note 
the pronounced curving of the streamlines as opposed to the rather straight, streaky 
pattern indicated by the oil film. The difference between the two patterns may be 
due to the fact that the oil film has a thickness and may be protruding into the 
boundary layer beyond the plane at which the numerical results are presented. The 
streamline pattern will show less curving as we move out of the boundary layer 
(see the pattern in Fig. 9 at an average distance from the surface of the wing 
d s = 6 x 10~ 3 c). 

Figure 10 shows the location of the separation line and reattachment line 
relative to the wing outline. Particles that start between the reattachment line 
and the wing tip are convected toward the tip vortex; particles outside that region 
will not be trapped by the tip vortex. The tip vortex can be clearly visualized by 
the three-dimensional path of the particles shown in Fig. 11. These particles are 
generated near the leading edge of the wing at several spanwise locations. Particles 
that start at a spanwise location y < 0.4 c do not get entrapped in the tip vortex, 
but particles starting at y > 0.4 c are captured in the tip vortex. The three- 
dimensional-flow separation on the wing surface forms part of the roots of the tip 
vortex. Figures 12 show the tip vortex as detected by particles released close to the 
tip of the wing. These figures show the projection of the tip vortex line on an x-y, 
x-z, and y-z plane. The spiraling of the particles around the wing tip is clearly 
seen from Fig. 12a and Fig. 12c. Figure 12b shows that the vortex lifts off the 
wing, the asymptotic angle reached is a = 5° (the angle of attack). 

The three-dimensional flow field is usually cumbersome to analyze, and an 
effective way to analyze such a field is to project surface sections of the field onto 
appropriate two-dimensional planes. Three surfaces, r) — 15, £ = 20, and £ = 
27, were selected to study the tip vortex formation. These surfaces intersect the 
surface of the wing, as shown in Fig. 10; that is, ?/ = 15 is a surface-cut of the 
flow field that will show the flow pattern in a section more or less parallel to the 
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streamwise direction; the other two surface-cuts show the flow pattern in sections 
approximately normal to the streamwise direction. Figure 13 shows the projection 
of the velocity vector at rj — 15 on an x-z plane. The flow- separation point and 
reattachment point are indicated by the arrows. The particle paths on this surface 
are shown in Fig. 14. Since the flow is three-dimensional, a recirculation zone is not 
necessary to indicate reattachment. The thickness of the separated layer is about 
0.075 c and is of the same order as the thickness of the boundary layer at the root 
of the wing; that is, the separation zone is fairly thin. Figures 15 and 16 show the 
projection on a y-z plane of the velocity vectors computed on the surfaces £ — 20 
and f = 27. The presence of a recirculation zone is more pronounced in these 
surface-cuts. Figure 17 shows the streamlines of the velocity vectors of Fig. 15. 
The core of the recirculation zone is clearly distinguishable in this figure. If this 
recirculation zone is taken to be the root of the tip vortex, we can conclude that the 
tip vortex is attached to the wing in this plane. The combination of the observations 
from the vector plots (Figs. 13 and 15) and the three-dimensional particle paths 
(Fig. 11) indicates that the vortex line on the wing surface away from the tip is 
almost parallel to the leading edge. Figure 18 shows the streamlines of the velocity 
vectors of Fig. 16. At this location, the core of the vortex is much smaller, and 
the center of the vortex has shifted toward the wing tip and lifted from the wing 
surface. The direction of this vortex is downstream. At this Mach number, we then 
have the picture shown in Fig. 19; that is, a vortex line parallel to the leading edge 
interacting with the tip vortex. 


YL Conclusions and summary 

A computer code has been developed to solve the three-dimensional, compres- 
sible thin-layer equations in generalized coordinates. The code can accommodate 
up to 200,000 grid points. A grid-generation system that is modular has also 
been developed that can be used to generate one-dimensional grid stretching, two- 
dimensional body-fitted grids, and three-dimensional body-fitted grids. This system 
was used to compute the flow field around a low- aspect-ratio wing at M = 0.82. At 
this Mach number, an oblique shock is observed on the upper surface of the wing. 
Also, a tip vortex forms and lifts off the wing. The formation of the tip vortex was 
analyzed. We found that the flow around the leading edge contributes to the vortex 
formation around the tip, and that the circulation generated in the boundary layer 
will also contribute to the formation of the tip vortex. 
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Fig. 2 Grid distribution. 
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Fig. 3 Grid wrap. 
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Fig. 4 Isobar lines on the wing surface. 
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Fig. 5 Mach contours on a surface parallel to the wing surface. 
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Fig. 6 Cp profiles at y — 0.6c. Comparison with experimental data 
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Fig. 8 Particle paths on f = 2 surface. 
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Fig. 9 Particle paths on f = 8 surface. 
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Fig. 10 Flow separation and reattachment lines. 
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Fig. 11 ' Three-dimensional particle paths. 
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(c) y-z projection 
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Fig. 13 Velocity vectors projection on an x-z plane, (surface r\ = 15). 
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Fig. 14 Particle paths on rj = 15 surface. 
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Fig. 16 Velocity vectors projection on a y-z plane, (surface f = 27) 
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Fig. 19 Sketch of the two vortex lines interacting on 
the upper surface of the wing 
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